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Analytical Representation for Equations of State of Dense Matter 

C. Giingor, K. Y. Eksi 
Istanbul Technical University 

Abstract. We present an analytical unified representation for 22 equations of state 
(EoS) of dense matter in neutron stars. Such analytical representations can be use- 
ful for modeling neutron star structure in modified theories of gravity with high order 
derivatives. 



1. Introduction 



The different possible equations of state (EoS) of dense matter inside neutron stars 
(NSs) are often given as tabulated data which are fed into the solution of the hydro- 
static equilibrium equations via an interpolation technique. While this method works 
well in the case of general relat ivity (GR) for which the hydrostatic equ ilibrium is de- 
scribed by the TOV equations ( Tolmanl l l939l : lOppenheimer & Volkoff 1939), it leads 
to pro blems in modified theories of gravity, e.g. f(R) theories dde Felice & Tsujikawa 
20101) in which the equatio ns describing the hyd rostatic equilibrium contain second 
derivatives of P(p) (see e.g. lArapoglu et al.l (1201 ll) ) where P is the pressure and p is the 
density. In such cases the order of the polynomial used in interpolation technique may 
effect the results if it is not sufficiently high and the linear interpolation technique will 
certainly fail. It is then necessary to employ analytical representations of the EoS that 
are sufficiently many times differentiable. In this work, we provide an accurate unified 
analytical representation for 22 EoS. 



2. Analytical Representation of EoS 



A unif ied analytical representation for two EoS, FPS and SLY, were provided by lHaensel & Potekhin 



( 20041) (hereafter HP04) where the authors provide a relation between £ = log(P/dyn cm" 
and £ = log(p/gcm -3 ) with 18 free parameters. The unified representation we provide 
here is an extension of HP04 with 23 parameters 12 o f which are fixed to represe nt low 



-) 



here is an extension ot HFU4 witn 15 parameters 12 of which are fixed to represent low 
density regimes of BPS dBavm et al.ll 19711) and NV dNegele & Vautherin|[l973l) . The 
rest 1 1 free parameters are used to fit the different EoS at high density regimes and to 
match them. The function we use to represent EoS for NS is 



t, = &>w/o(«l(£ - en)) + Ma 2 (cn ~ £>Khi g h 



(1) 



where fo(x) = 1/(1 + exp x) is the matching function (also used in HP04). 
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Figure 1. The EoS AP4 (+) and its ana- 
lytical representation (red line). The purple 
dotted line and turquoise dash-dotted line de- 
scribe low and high density regimes, respec- 
tively. The small figure shows the adiabatic 
index F = d£/cUj versus 



Radius (km) 

Figure 2. M-R relation 
for EoS, AP4 and MPA1. 
The red line is obtained 
by the analytical represen- 
tation of the EoSs as pre- 
sented in this work while 
the blue circles are ob- 
tained by the tabulated 
EoS. 



Here 

£ow = [c\ + c 2 (t - c 3 ) C4 ]f {c 5 {£ - c 6 )) + (c 7 + c^)f {c 9 (c w - £>) (2) 

and 

4igh = («3 + a^)fo («5 (^6 - 0) + (an + a%£ + a 9 ^ z )f (a l0 (a n - 0) 0) 

describe the low and high density regimes, respectively. We provide the values of the 
fit parameters c,- and a, for £ > 5 in Table Q] and Tabled 

In FigureQ] we show our results for fitting the EoS data of AP4 with Equation £T|). 
In Figure |2j we compare the M-R relations obtained by the analytical representation 
with that obtained by feeding the EoS data via interpolation technique. The maximum 
relative error is ~ 0.05% near the maximum mass. 



3. Discussion and Conclusion 



The unified analytical expression presented in this work is an accurate representation of 
many EoS suggested for NSs. The analytical representation is preferable for solving the 
structure of NSs in modified theories of gravity where hydrostatic equilibrium equations 
are of 4th order. In such cases the usual interpolation technique will fail because the 
high order derivatives may not be continuous if the order of the polynomial used for 
interpolation is not sufficiently high. 
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Table 1 . Parameters of the fit for low regimes 



Cl 


10.6557 


C2 


3.7863 


C 3 


0.8124 


C 4 


0.6823 


C5 


3.5279 


C6 


11.8100 


C 7 


12.0584 


c 8 


1.4663 


c 9 


3.4952 


ClO 


11.8007 


Cll 


14.4114 


C12 


14.4081 



Table 2. Parameters of the fit 



EoS 


API 


AP2 


AP3 


AP4 


engvik 


gmlnph 


gm2nph 


gm3nph 


mpal 


msOO 


ms2 


ai 


4.3290 


4.3290 


6.3293 


4.3290 


5.0487 


11.4832 


l A fin a 1 

10.9801 


11.5163 


5.2934 


16.6491 


14.0084 


a 2 


4.3622 


4.3622 


6.3467 


4.3622 


5.0838 


11.4006 


10.9233 


11.4210 


5.3319 


16.2520 


13.8422 


a 3 


138.1760 


184.1600 


174.2030 


9.1131 


-100.7550 


-100.7660 


-100.9130 


-100.8670 


87.7901 


8.0809 


16.5970 


a 4 


-10.1093 


-13.4285 


-15.0960 


-0.4751 


9.2091 


11.8432 


15.6165 


16.0001 


-5.8466 


-0.5036 


-1.0943 


a 5 
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J.V / Ul 
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14.0120 


13.8551 


13.8962 


14.8800 


13.9137 


12.7088 


11.9067 


11.7186 


15.0804 


14.7513 


14.8169 


a 7 


-411.1380 


-372.6920 


-377.4200 


21.3141 


462.5180 


462.0230 


461.2360 


461.2870 


428.4130 


-39.9099 


-56.3794 


a 8 


48.0721 


40.4186 


44.3688 


0.1023 


-48.4401 


-51.9392 


-55.7395 


-55.2969 


-57.6403 


8.0314 


9.6159 


a 9 


-1.1630 


-0.8061 


-0.8911 


0.0495 


1.1892 


1.2438 


1.2502 


1.2024 


2.0957 


-0.2022 


-0.2332 


aio 


4.7514 


4.7456 


4.5062 


4.9401 


5.0896 


5.0896 


5.0897 


5.0897 


5.0588 


5.0102 


-3.8369 


an 


10.7234 


10.7503 


13.6882 


10.2957 


10.2277 


10.2277 


10.2276 


10.2276 


10.2727 


10.3408 


23.1860 




msl506 


pal2 


pclnphq 


wffl 


wff2 


wff3 


wff4 


schafl 


schaf2 


prakdat 


ps 


ai 


15.4744 


11.3079 


11.3180 


2.8144 


3.9694 


5.2934 


5.0395 


11.9204 


13.3358 


12.3641 


12.3769 


a 2 


15.1705 


11.1963 


11.2234 


2.8539 


4.0017 


5.2912 


5.0512 


11.7992 


13.1285 


12.2125 


12.0277 


a 3 


2.8727 


-0.3393 


-9.2769 


-6.6340 


-6.1548 


-5.8771 


-5.9160 


-11.7736 


-6.0249 


-9.1898 


-4.5181 


a 4 


-0.2014 


0.0149 


0.6194 


9.3448 


20.5349 


29.4964 


29.5087 


3.1670 


2.7698 


1.0249 


2.6750 


a 5 


4.2468 


5.0914 


7.9411 


0.8728 


0.5724 


0.2267 


0.1944 


-12.0189 


-13.5535 


2.5784 


-11.4273 


a 6 


15.1919 


14.5066 


15.2677 


15.2163 


15.3019 


17.8721 


21.5182 


14.6274 


14.7235 


13.6433 


14.5886 


a 7 


-83.9695 


-68.4757 


-194.2320 


-379.9680 


-377.0480 


-380.2820 


-381.1690 


15.3263 


-398.4300 


44.5552 


-84.9113 


a 8 


13.8824 


11.8013 


28.9879 


77.0628 


85.8006 


65.3388 


61.2652 


0.5165 


55.1217 


-2.3817 


13.9429 


a 9 


-0.3959 


-0.3266 


-0.9137 


-3.5619 


-4.5011 


-3.1750 


-2.6609 


0.0521 


-1.7487 


0.0902 


-0.3964 


aio 


-4.3087 


-4.3089 


-4.3181 


-4.3181 


-4.3181 


-4.3181 


-4.3181 


12.1722 


13.6495 


4.8729 


11.5048 


an 


23.3603 


23.3604 


23.3626 


23.3626 


23.3626 


23.3626 


23.3626 


14.6270 


14.7214 


10.3969 


14.5874 



